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We consider a graphene sheet in the presence of a strong perpendicular magnetic field with a single 
short-range 5-impurity situated at one of the carbon sites. We study the neutral inter-Landau 
level collective excitations, magnetoplasmons, which become localized on the impurity. Some of 
these excitations involve a pseudospin flip (intervalley transitions), since the impurity can scatter 
electrons between the two valleys. We propose a classification of states of the excitations in graphene 
and introduce the appropriate quantum numbers. The energies and optical strengths of collective 
excitations are calculated for a range of integer filling factors and impurity strengths. We establish 
a set of symmetries matching the energies and absorption strengths of collective excitations for 
different sublattice locations of the impurity, filling factors, circular light polarizations and signs of 
the impurity potential. 
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I. INTRODUCTION 

At first glance, the electronic structure of graphene 
is simple: due to the sp 2 bonding of the single-atom- 
thick two-dimensional (2D) carbon layer, it is easily mod- 
eled by a non-interacting tight-binding approach*^ For 
the 2D honeycomb lattice, this results at low energies in 
two non-equivalent Dirac cones at the K and K' points. 
And indeed much of the tremendous progress 3 ^ in un- 
derstanding graphene's electronic properties^ since its 
isolation^ in 2004 has been due to the exploitation of 
this fact. However, electron-electron (e-e) interactions^ 
also play an important role in graphene, particularly in 
the high magnetic field regime^— As well as the trans- 
port properties, e-e interactions affect the optical excita- 
tions of the system. Previous studies have calculated 
the dispersion relation for particle-hole excitations for 
the case of clean graphene with integer Landau level 
(LL) filling ii3ri& This has been done for the bilayer sys- 
tem tooj 16 ' 17 Furthermore, charged collective excitations 
have been predicted to exist as discrete states outside the 
continuum ji^ 

In this work, we investigate the effect of short-ranged 
disorder on the collective excitations (CEs) of the mono- 
layer system in the presence of a strong perpendicular 
magnetic field. Specifically, we calculate the energies and 
optical properties of CEs that become localized on a 6- 
function scatterer situated at one of the graphene lattice 
sites. We shall refer to this scatterer as an impurity, but 
it can represent a vacancy as well . 19 ' 20 We explore how 
the bound states arc influenced by the filling factor v, 
the sublattice containing the impurity, the light polar- 
ization and the impurity strength. The results are no- 
tably different from those for a long-range Coulomb im- 



purity potential, studied previously^! In particular, the 
choice of a sublattice position for the impurity breaks 
the equivalence between K and K' points and thus the 
SU(4) symmetry. Furthermore, the short-range character 
of the impurity allows for significant inter- valley scatter- 
ing. Nevertheless, we will show that due to particle-hole 
symmetry at the Dirac point, the resulting optical excita- 
tions obey another form of sublattice symmetry between 
attractive and repulsive impurities. Our results show a 
symmetry structure which allows us to distinguish be- 
tween long- and short-ranged impurities in the optical 
excitations of graphene. 



II. THEORETICAL APPROACH 

A. Single particle problem 

Let us consider a single electron in graphene in a per- 
pendicular magnetic field with a (5-function impurity lo- 
cated at the origin. For the case when the origin is at a 
H = A,B sublattice site, the Hamiltonian has the form 
H^ = Ho + Vfj,, where Hq is the Hamiltonian for a free 
electron in graphene in a uniform magnetic field 3 
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Here vp is the Fermi velocity and n± = n 
II = p + -A, the kinematic momentum operator. Note 
that the (AK, BK, AK' , BK') ordering is used and that 
Hq is diagonal with respect to the valley index. The 
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FIG. 1: (a) Real space graphene lattice with nearest neighbor 
vectors defined for the two cases of having anAorB sublattice 
site at the origin, (b) Reciprocal lattice with the first Brillouin 
zone indicated by shading and definitions of K, K' points. 



contribution to the Hamiltonian from the short-range 
V(r) = VS (r) impurity is 



V a = V6(t) 
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for an impurity located on an A site at the origin and 

u ; = ivnr)i ;; /, ;; ;; i. (3) 



for an impurity located on a B site at the origin. Here 
a = e 27 ™/ 3 and V — y/SWa 2 /2, where W is the onsite 
energy associated with the impurity, see Eq. (IA1|) : a is 
the distance between atoms on the same sublattice, see 
Fig. Q] The off-diagonal terms in Eqs. © and © de- 
scribe scattering between the valleys; such terms can be 
neglected for long-range potentials. Note that in both 
cases n = A,B the origin is chosen at the impurity. The 
derivation of Eqs. ©-© is given in Appendix [S] see 
also Ref. [22|. The eigenstates of the Hamiltonian Hq in 
the symmetric gauge A = ^B x r for an electron in, e.g., 
the K valley (pseudospin r = ft), are 
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Here n is the integer LL number, 4>\ n \ m (r) is a 2D Elec- 
tron Gas (2DEG) wavefunction with oscillator quantum 
number m = 0, 1, . . ., a n — 22( <5 >»>°~ 1 ) ; S n = sign(n) (with 
<So = 0) and % s is the spin part for the two spin projec- 
tions s = f,4- The corresponding wavefunction in the 
K' valley (r = 4) is obtained by reversing the order of 
the spinor components. There is the Landau degeneracy 
in quantum number m, so that defining the composite 
index Af = {nsr}, the single particle energy is 



(-M = S n huj c 
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Here hw c — \/2hvp /£b is the cyclotron energy in 
graphene, £b = \JhcjeB is the magnetic length and 
fvjj s , hcu v the phenomenological spin and valley splittings, 
respectively. We assume hio s > Hu) v and that these split- 
tings are small in comparison with the cyclotron energy, 
hjjj s ,hui v <C hu) c ; we set them equal to zero for the pur- 
poses of numerical calculations. Note that the wavefunc- 
tions \l/_A/ m (r) are not eigenfunctions of the orbital angu- 
lar momentum projection operator l z = r x p. Instead, 
the generalized angular momentum 2 ^. j z = l z + \(T Z ® I is 
conserved in graphene. Here a z is a Pauli matrix acting 
on the isospin (in the sublattice index space /i = A, B) 
while the unit matrix I acts on the pseudospin (in the 
valley K,K' space r = ft, if). The generalized angular 
momentum operator gives j z ^Nm{^) = jz^Mm( r )i with 
the half-integer eigenvalue j z = l z — h and the integer 
orbital part l z — \n\ — m. 

The calculation of impurity matrix elements involving 
Vn is straightforward, since 
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where / <G {|n|, \n\ — 1} and k £ {\n'\, \n'\ — 1}, with ex- 
act values set by pseudospins (valley indices) r and r'. 
In addition, 4> nm (0) ~ S ntTn , imposing the selection rule 
determining which states are affected by the (5-impurity. 
Indeed, only the s-orbitals with l z = \n\ — m = have 
non- vanishing probability amplitudes at the origin, where 
the <5-impurity is located. Note that this selection rule 
allows mixing of states with different generalized angu- 
lar momentum projections j z . This is in contrast to a 
Coulomb potential, where j z is strictly conserved. In 
Appendix [B] we dicuss when higher order corrections to 
the energies ejv due to the impurity interaction are small 
and may be ignored. 



B. Collective excitations 

We consider excitations in which an electron is pro- 
moted from one of the uppermost filled LLs, A/2 in the 
Dirac sea, to an empty higher lying LL A/i, leaving be- 
hind a hole. The creation operator for such an excitation 

is 
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where the hole representation, cj^ m — > d\f m and cL- m — » 
dj\f m , is used for all filled levels. The excitation oper- 
ator QJyy ^ j acts on the ground state, denoted by \v) 
for a system with filling factor v. The expansion co- 
efficients satisfy ^^^(toi, m 2 ) ~ 5 Jjti | ni |_ mi _| na | +rna . 
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FIG. 2: (Color online) Particle-hole excitations for filling fac- 
tor v — — 1. The thin arrows indicate two possible spin states 
s = f, 4- and the thick arrows two possible pseudospin r = "ft, JJ. 
(valleys K, K') states. The dashed lines connect excitations 
mixed by the direct Coulomb interaction, see Fig. [3ji. The 
dashed box contains four excitations resonantly mixed by the 
exchange Coulomb interaction, see Fig. [3)3. The six leftmost 
(red-dotted) excitations are those which are mixed when there 
is an impurity on the A sublattice and are optically-active in 
the left circular polarization a + of light. 



The quantum number J z is the total generalized angu- 
lar momentum projection given by summing over the 
j z values for individual particles i, j Z i = l z % — \, with 
hi = \n-i\ — mi- In second quantized form, 
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the neutral particle-hole excitations J z — l z \ — l z % — 
\n\\ — m,\— \ri2\+m2 only contains an orbital part, in con- 
trast to the single particle states and charged collective 



excitations^ in graphene. 

The total Hamiltonian, H^ = He^+Hh^+H-mt, includ- 
ing the free energies, the interaction with the impurity on 
the fi = A,B site, and the electron-hole (e-h) interaction 

is 

H H = E [^MxM^rmm'^Mi + Hwimi J C M[ m ' 1 C M 1 m 1 

Mi Ml 
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The ew are the single particle energies, which are renor- 
malized due to the exchange interaction with other elec- 
trons in the Dirac sea (see Appendix [C|l . The last term 
gives the dynamical part of the e-h Coulomb interaction, 
following from the pairwise Coulomb potential of e-e in- 



teractions in graphene, f7 (|ri — r 2 | 



Here 



£|ri-r 2 | ' 

e is an effective dielectric constant, which depends for 
graphene on its environment. The dynamical e-h inter- 
action is made up of the e-h direct attraction and the 
e-h exchange repulsion H lrA = H® h + H* h with the total 
vertex given by 
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see Fig.[3J The matrix element U is defined in the electron 
representation by 
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wfiSS = S. lA 8. aA /dri|dr 2 ^^(n) ® ^ m ^ 2 )U (|n - r 2 |) * nirtmi (n) ® $„ 2T2m2 (r 2 ) , (11) 

I 



where ® denotes the direct (Kronecker) product. Note C. Symmetry properties of the Hamiltonian 

that IA may be expressed in terms of the correspond- 
ing matrix elements for the 2DEG2i and conserves spin, 
pseudopsin (no intervalley scattering), and the general- 
ized angular momentum: 
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It also possesses the SU(4) symmetry,—^ 
below. 
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as described 



In the absence of an impurity (V = 0), the Hamiltonian 
H = Hq + Hi n t has several symmetries. Firstly, there 
is the generalized axial symmetr y 21 1 23 [J Z ,H] = 0. This 
symmetry is broken in the presence of a <5-impurity, which 
scatters particles between the valleys. Specifically for 





(a) 






niSiTimi 

> — 






n'jSiTiTTl'j 
> 


■*— 

n 2 S2T 2 m 2 






< 

n'. 2 S2T 2 m'. 2 




(b) 


Heh 




msirimi 






n\s' l T' 1 m! 1 




ri2S\T\m 2 



n'nS't TiTTl'.-, 



FIG. 3: Direct (a) and exchange (b) electron-hole interaction 
vertices, see Eq. (|10[) . 
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withj^j— j z i = |n^| — mj — |ni|+mi. This results in excita- 
tions with different J z quantum numbers being mixed by 
the impurity, which we shall discuss in Sec. Ill El There- 
fore, J z ceases to be an exact quantum number. We will 
still use it for classification of states to indicate their ori- 
gin and predominant character. 

Another kind of symmetry present in clean graphene 
with long-range Coulomb interactions is the SU(4) sym- 
metry, due to the equivalence of the four possible 
pseudospin-spin states, UTi-UTifNU-O-l- The generators of 
the SU(4) group may be expressed in our context in sec- 



ond quantized form as C jk = J2n c N 3 c Nk - Y^n <F Nk d Nj , 
where N = {nm} is a set of orbital quantum numbers 
and j, k = 1rt,JJ-t,1U,-U-l arc the f° ur SU(4) flavors. In 
the absence of an impurity and for a fully filled LL, 
[H, Cjk] — 0. Upon introducing an impurity, the com- 
mutation relation becomes 



[H„c Jk \ = J2 E( 1+, W{) 



(14) 
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For an impurity which does not change the flavor of a 
particle and has scattering matrix elements which are 



N'f _ 



flavor independent, V M7V / = SkjV^j^ , one can see from 
Eq. (TT4")) that the SU(4) symmetry is still preserved, 
[H/j,, Cjk] = 0. This is the case for, e.g., a Coulomb impu- 
rity and explains the degeneracies of the corresponding 
states^ In contrast, the <5-impurity considered here may 
scatter between the valleys, thus changing a particle's 
flavor and breaking the SU(4) symmetry. 



D. Optical selection rules 

We work in the electric dipole approximation, ignor- 
ing the magnetic field component of the electromagnetic 
wave. The Hamiltonian describing the electron-photon 
interaction for an incoming circularly a ± polarized beam 
of light is 



8H ± = l ^lL e -^ 



o± 

CTzr 



(15) 



where T represents the electric field strength, u is the an- 
gular frequency and a± — a x ± io~ v are the Pauli matrices 
acting on the isospin (sublattice A, B) components. In 
the dipole approximation, the linear momentum trans- 
ferred to the electron by a photon is negligible; accord- 
ingly, no intervalley transitions can be induced. Besides, 
the electric field is not coupled directly to the electron 
spin, so that 5'H± conserves spin, i.e., no spin flips occur 
in electric dipole optical transitions. In addition, one can 
show that 



(Af'm'\Sn±\N'm) ~<5 



KlTi,M c 



where \Afm) = cj^ m |0). The single particle selection rules 
are thus \n'\ — \n\ = ±1, to = to', t = t' and s = s', 
for the a polarizations, where the unprimed quantum 
numbers describe the initial state of the electron and the 
primed quantum numbers its final state, see Refs. i21,26l , 
[271 and references therein. 

In the second quantization form and in the e-h repre- 
sentation, the interaction with the circularly a polarized 
light is 



SH ± = J2 (n'\SH±]n)J2' 
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Here we made use of the fact that the matrix element 
{N'm\5i-L±\Nm) = 5 s ^ s 5 T ^ T 8 m ^ m {n'\5'H ± \n) is diagonal 
in s, t, to and, besides, it does not depend on these quan- 
tum numbers. The operator emerging in Eq. (|16l) . 
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is proportional to the composite boson creation operator 
Q\f N k=o describing a magnetoexciton of zero magnetic 
momentum k = with electron (hole) in the N' (M) LL 
in graphene. Nq — S/2nl 2 B is the LL degeneracy with 
S the area of the graphene sheet and the operator of 



magnetic translations is k = X^C^-i — — i"i x B). More- 
over, the electron-photon interaction is SU(4)-symmetric, 
[Sy.±,Cjk] = 0; it does not change the flavor of an elec- 
tron and has matrix elements which are flavor indepen- 
dent. This leads to the coupling of the electric dipole 
photon only to the flavorless boson in the given LLs2i 

Ql'n = £ s ,t Ql' S msrn=o- In the presence of an impu- 
rity, the magnetic momentum k is no longer conserved, 
[Hfj,,k] ^ 0. In addition, a short-range impurity induces 
the intervalley scattering and breaks the SU(4) symme- 
try, further relaxing the selection rules. 

Let d = {N\N 2 J z \8'H±\i>) denote the optical dipole 
transition matrix element from the ground state \v) with 
filling factor v, to the final state with one excitation 
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The transition matrix element squared |d| 2 is propor- 
tional to the absorption intensity or oscillator strength 
for that state. It can be shown that 
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so the optical selection rule for a CE is J z — ±1. Besides, 
only the excitations with the total spin and pseudospin 
projections S z = s zl - s z2 = 0, T z = t z1 - t z2 = 
are optically active. The latter selection rule is strict. 
For equal spin (pseudospin) filling of a given n th LL, the 
CE states can additionally be classified by a total spin 
S = 0, 1 (pseudospin T = 0, 1). Among those states only 
the spin S = and pseudospin T = singlets have non- 
zero projections onto the state Q n / n , which is directly 
coupled to photons. Therefore, only the singlets S = 0, 
T = are optically active, while all the triplet S = 1 and 
T = 1 states are optically dark. 



E. Mixing of excitations 

The long range nature of the Coulomb potential means 
it cannot provide a large enough change in momentum 
|Ak| ~ |K — K'| ~ 1/a to scatter between the valleys. 
It also only mixes excitations with the same angular mo- 
mentum projection J z = |m| — mi — \n2\ + Tn<i. In ad- 
dition, Eqs. (fTU)) and (JT2J) restrict the possible spin and 
pseudospin states of transitions mixed by the Coulomb 
interaction. Generally, two excitations, \MiN 2 J z ) an d 
l-A/pV^ J z ), are mixed by (i) the direct interaction if 
s\ = s 1 , n = r 1; s 2 = s 2 and t 2 = t 2 (see Fig. [SJi) and 
(ii) the exchange interaction if si = s 2 , T\ = t 2 , s 1 = s 2 
and Tj = r 2 (see Fig. [3b) . Figure [2] shows the possible 
excitations for filling factor v = — 1. The four excita- 
tions with no spin or pseudospin flips contained within 
the dashed box are mixed by the exchange interaction. 
Pairs of excitations connected by a dashed line are mixed 
by the direct interaction. The remaining excitations are 
unmixed by the e-h Coulomb interparticle interactions. 

In our calculations we assume all LLs with n < are 
filled, all LLs with n > are empty and that the sublevels 



of the zeroth LL become successively completely filled 
(v = —1,0, 1,2). We have seen that infinitely many ex- 
citations \N\N 2 J Z ) with the same J z and particular spin 
and pseudospin quantum numbers are mixed by the e-h 
Coulomb interaction. However, we may truncate the ba- 
sis so as to obtain a tractable finite-size matrix represen- 
tation of Hamiltonian H^, Eq. ([9]). To this end we only 
consider excitations, which are either in resonance (have 
same energies) or are very close to resonance. In particu- 
lar, we are interested in the lowest energy ~ huj c inter-LL 
excitations for the chosen filling factors. Therefore, we 
need to consider the mixing of — > 1 and — 1 — > ex- 
citations. Neglecting the mixing of non-resonant excita- 
tions amounts to ignoring terms with amplitudes, which 
depend on the "fine structure constant" for graphcnc, 
a = e 2 jhvpe. 

But even with this simplification, the Hamiltonian ma- 
trix remains infinite because of the macroscopic degen- 
eracy of each LL in the oscillator quantum number m. 
However, as long as we are interested in the states lo- 
calized on the impurity, we may truncate the summation 
over m in Eq. ([7]) at a finite sufficiently large value m max . 
This is justified because the distance from the origin (im- 
purity) of single-particle orbitals in LLs increases with m 
as 



(r 2 ) = (2\n\ + 2m+l + 5 nfi )£ 2 B . 



(20) 



As a result, for e-h states localized on the impurity we 
have a fast convergence with mi 21 ' 28 In our numerical 
calculations we choose m max = 30; the total number of 
states involved ranges from 180 to 240, depending on v, 
the sublattice of the impurity and which light polariza- 
tion the excitation is bright in. 

We are primarily interested in optically active exci- 
tations, although these will always be mixed by the 5- 
impurity potential to some which are nominally dark. 
This leads to the redistribution of the oscillator strength 
so that a number of states instead of being strictly dark 
become "faint" . In the previous Sec. Ill PL we saw that 
bright states had no spin or pseudospin flips and J z = ±1 
for the a polarizations. The four excitations with no 
spin or pseudospin flips are always mixed by the exchange 
e-h Coulomb interaction, although how many of them are 
bright depends on the filling factor and light polarization. 
Two or four additional excitations with pseudospin flips 
will be mixed to these by the (5-impurity; generally, for a 
given J z , the admixtures will have J z and J z ± 1 values. 
Thus the full CE comprises of six or eight excitations 
with different spin/pseudospin characters and J z values. 
As an example, Fig. [2] indicates the six (red) excitations 
which should be taken into account for v = — 1, the a + 
light polarization and an impurity on the A sublattice. 
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FIG. 4: (Color online) Dimensionless energies £ of localized 
CEs vs dimensionless impurity coupling constant g [Eqs. (I22|l . 
(|23|l ]. for the (A, a + , v — — 1) system (squares) and the con- 
jugate system (B, a~ , v — 1) system (circles). Notice the 
symmetry of the spectra, Eq. (|26[) . The shaded region corre- 
sponds to the magnetoplasmon band with extended states. 



III. RESULTS AND DISCUSSION 

A. Energies and optical strengths 

A representative selection of results are displayed in 
Figs. HH3 In the absence of an impurity / 3 ! 14 ! 25 all CEs 
have extended wavefunctions and fill a band (indicated 
by the shaded regions) of width determined by the char- 
acteristic Coulomb energy 
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B. 
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More precisely, the band has width Eq for odd filling fac- 
tors v = —1,1 and width 0.75Eq for even filling factors 
v — 0, 2. The difference occurs because excitations which 
are mixed by the direct interaction, such as the pair in 
Fig. [2] with J z = 0, partake in the CE for odd filling 
factors, but not for even ones (see Appendix [C]l . Upon 
introducing an impurity, discrete localized CEs with en- 
ergies E appear both below and above the band of con- 
tinuum states. We plot the dimensionless energy £ as a 
function of the dimensionless impurity coupling constant 
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(22) 
(23) 



Thus the impurity strength W and magnetic field B may 
be simultaneously tuned so that the quantity g is left un- 
changed. Note that in this case the dimensionless ener- 



gies £ remain unchanged, indicating a scaling property 
of our theory. 

Note that all the plotted energies have been renormal- 
ized: The lower continuum edge corresponds to the ob- 
servable lowest cyclotron mode, h£j c = ftio c + 5huj c , which 
has been set to the zero energy reference level. Here 
0.7-Eo is the e-e correction to the 



Shur. = Aitl + T m 



"Oft 



on 



bare cyclotron energy, with A H ~ 1A3Eq the total self 
energy correction and r o H = — 0.75E the "excitonic" or 
vertex correction (see Appendix [Cl and Ref. 21). The size 
of square data points in Figs. 2H7]is proportional to the 
optical dipole transition matrix element squared |d| 2 , see 
Sec. Ill Dl above. Notice that the majority of states are 
optically active, although some (marked as "faint") are 
considerably weaker than others, by at least three orders 
of magnitude. Some branches become brighter as they 
approach the band, whereas for others, their oscillator 
strength decreases. The physical reason for this behavior 
is at present unclear. 

The dark states shown in Fig. 0] have a special charac- 
ter. Out of the six possible types of excitation in our basis 
(see Fig. [2]), only the two excitations, |04-IJ- — 144 J z = 1) 
and |0t-ll- -It! Jz = 1) contribute; see Eq. (JTSJ for 
the definition of |A/]7v2</z)- Having no amplitude on 
|lt-D- Ot-IJ- J z = 1), the only excitation out of the six 
which is optically active in o~ + , explains why they are 
dark. Explicitly they are created by the operator 



Dt 



>(m)( 



r 1 d 



-lt-U-2+m 



r ! d 



•1442+m 



(24) 
with certain amplitudes A(m) rapidly decreasing with m 
[distance from impurity, Eq. (I2U)) ]. Although generally 
the total spin quantum number is not well defined for 
CEs, excitations (|24[) have a higher symmetry and are in 
fact spin triplet states S = 1,S Z = 0. Besides, the type 
of excitations is unusual in the following sense: All other 
excitations shown in Fig.|4]have contributions from all six 
of the excitations which can be mixed. We note that the 
excitations with pseudospin flips are in general strongly 
mixed by the (5-potential to those with no pseudospin 
flips. 

The results indicate critical g values for the formation 
of localized excitations, although within our approach it 
is difficult to resolve what may be a barely bound low- 
energy state from the lower continuum edge. The crit- 
ical value for which bound states appear varies slightly 
depending on the filling factor v : light polarization a 
and the sublattice position A, B of the impurity. For 
example, for a system with filling factor v = — 1 illumi- 
nated by er + circularly polarized light with an impurity 
on an A sublattice [such a system will henceforth be de- 
noted by (A, cr + , v = — 1)], there are no bound states 
for -2.3 < g < 7.4. If we then take B = 15 T and e = 5 
for example, this corresponds to an impurity strength 
of approximately -100 eV < W < 340 eV. Notice that 
short range potentials with rather large amplitudes cor- 
respond to vacancies^ It takes a very strong (5-function 



impurity potential to localize excitations, since the impu- 
rity only couples to basis states where one of the particles 
has m = 0, 1. Strictly speaking, such energies are larger 
than the width of the it band, where the electrons can be 
reasonably treated as massless Dirac fermions. For such 
high values of the impurity potential, our results have 
only a qualitative nature. 

In all cases, larger magnitudes of g are required to push 
the states above the band than to pull them below. This 
is because the states below the band form when either the 
electron or hole is nearer to the impurity and attracted 
to it; in this case the additional e-h attraction also low- 
ers the energy. States above the band form when one 
particle is held nearer to the impurity due to magnetic 
confinement, but is repelled by it, whilst the other parti- 
cle is further away. In this case the e-h attraction works 
against this, trying to lower the energy, so that a larger 
impurity strength is required to overcome this. 

The bound states move further from the band as the 
impurity coupling constant g increases, as expected. An 
interesting question is what happens to the bound states 
when they approach the band. One possibility is they 
cease to exist, merging with the two-particle e-h con- 
tinuum. The alternative is that some continue to exist 
within the band as quasibound states (resonances) . As a 
general rule, the latter states have high probability am- 
plitudes on the impurity and long-range oscillating tails, 
which make them non-normalizable^ The existence of 
resonances is facilitated by the confining effect of the 
magnetic field. Although we detect possible signatures of 
such resonant states, our method is not accurate enough 
to claim their existence. Interestingly, in all cases, the 
number of branches below the band for positive (nega- 
tive) g equals that of branches above the band for nega- 
tive (positive) g. 

In most cases for large enough \g\ values, bound states 
may be found both above and below the band for both 
g > and g < 0. One exception (see Fig. [4]) is the (A, 
o~ + , v = — 1) system. In this case the impurity selec- 
tion rule forbids the hole to interact with the impurity, 
so that all diagonal matrix elements have the same sign 
as g. As a result, the states above the band are only seen 
for large enough positive g and states below the band 
are only for large enough negative g. This behavior is 
mimicked by the "conjugate" -symmetric (see Sec. IIIIBI 
below) system (B, <r~ , v = 1) except that the electron 
is forbidden to interact with the impurity. As a result, 
states above the band only appear for large enough nega- 
tive values g < and states below the band only for large 
enough positive g > values. The other exception is the 
pair of "sublattice" symmetric systems (A, a~ ,v = 2) 
and (B,a~ ,v = 2), where in both cases only the hole 
interacts with the impurity so that states above (below) 
the band only exist for large enough negative (positive) 
values of g. These are in turn "conjugate" -symmetric 
to the systems (B 7 a + ,u = —2) and (A,a + ,v = —2) re- 
spectively. The number of branches due to the electron- 
impurity and the number of branches due to the hole- 
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FIG. 5: (Color online) Energies of localized CEs vs impurity 
strength as in Fig. U but for the symmetric systems (A, a + , 
v — 0) and (B, a + , v — 0). States marked as "faint" have 
values of |d| 2 at least three orders of magnitude smaller than 
the brightest states. 
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FIG. 6: (Color online) Energies of localized CEs vs impurity 
strength as in Fig. \S\ but for the (B, a + , v — — 1) system. 



impurity interaction can be correctly predicted by study- 
ing the transitions involved and the impurity selection 
rules. 



B. Symmetries of the spectra 

The above results indicate that the CEs possess a high 
level of symmetry. Moving the impurity from one sub- 
lattice to the other, A ■<-> B, yields the same energies 
and oscillator strengths for filling factors v — 0, 2 (see 
Fig. [5]). However, the results are qualitatively different 
for filling factors v = ±1 (compare Figs. 0] and [5]). This 



is because for odd filling factors such as v = ±1, the val- 
leys and thus sublattices are inequivalent. This follows 
from the phenomenologically introduced valley splitting 
energy Hu) v t z [see Eq. ||5J)], which leads to unequal oc- 
cupancies of the two valleys. However, for even filling 
factors v = 0, 2, . . ., the occupancies are the same and ef- 
fectively there is no valley splitting. For the cases when 
the sublattices may be considered equivalent, the A ^ B 
correspondence of the energies and oscillator strengths is 
understood formally by observing that the Hamiltonians 
Ha and Hb are connected by a unitary transformation 
Hb = UHaU 1 ' ■ This transformation interchanges both 
valley and sublattice indices and is given by 



U \ n _™* n n > ( 25 ) 












t) 








a 





—a* 








a* 








0/ 



with a = e 27 "/ 3 . 

Upon moving the impurity between the sublattices 
i f> B, symmetry may be restored also for odd fill- 
ing factors. This requires the valleys to be interchanged, 
which can be achieved by reflecting about n = 0. This 
entails flipping the spin and pseudospin quantum num- 
bers and transforming electrons into holes and vice versa 
(n — > — n), which means that v ■<-> —v and the sign of the 
impurity potential and also the light polarization should 
be changed. Note that the LL structure in graphene, as 
well as its upper-lower cone dispersion in the absence of 
a magnetic field, is electron-hole symmetric. As a result, 
the energies and oscillator strengths of CEs induced by 
circularly polarized light in a system with filling factor v 
and an impurity on the A sublattice of strength V, cor- 
respond exactly to those of excitations induced by light 
circularly polarized in the opposite direction in a system 
with filling factor — v and an impurity on the B sublat- 
tice of strength — V. This symmetry may be expressed 
symbolically as 



A, v, V, cr d 



B.-z/.-F,^ 



(26) 



An example is shown in Fig. 2) The above symme- 
try can be seen formally as follows. Switching the 
filling factor v — > —v changes a pair excitation 33 



\Af 1 m u Af 2 m 2 ; v) = c\ /im d\ f , m _\u) according to 



|A/imi,A/2m 2 ;^) ->• \M2m2,Mimi', — v) 



(27) 



where M = {n,s,r} and the conjugate M = {— n, s, f} 
with s, f representing flipped spin and pseudospin. From 
this it follows that J z — > — J z . From the e-h symmetry, 

(-v; N! 2 m 2 ,N[m' l \H in tW2m2i Nina; -v) = 
= (iy;M[m' 1 ,^m' 2 \H int \M 1 m 1 ,M 2 m 2 ;iy)* . (28) 

Hence the two matrices have the same energy eigenvalues 
and conjugate eigenvectors leading to the same values of 
the dipole transition matrix elements \d\ 2 . 
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FIG. 7: (Color online) Energies of localized CEs vs impurity 
strength as in Fig. [S] but for the (A, a~ , v = —1) system. 



IV. CONCLUSIONS 

We have considered binding of neutral collective exci- 
tations on a short-range impurity in graphene in strong 
magnetic fields. The considered excitations are formed 
when an electron is promoted to one of the empty Lan- 
dau levels leaving behind a hole in a lower filled level. 
This kind of problem requires treatment of interparti- 
cle (electron-hole) Coulomb, electron-impurity, and hole- 
impurity interactions all on equal footing. The scheme 
developed in this paper treats the impurity in the con- 
tinuous approximation, can be applied to arbitrary inte- 
ger LL filling factors v and gives the energies and opti- 
cal strengths of excitations localized on the impurity in 
graphene. The concrete results we have presented and 
discussed here are for the n = Landau level filling fac- 
tors v = -I, 0,1,2. 

A separate study^S using the same approach has shown 
that the <5-impurity localized states do not seem to evolve 
smoothly as a function of filling factor v in contrast to 
the case of the Coulomb impurity. 2 ^ This is perhaps be- 
cause for the Coulomb impurity there are always four 
transitions with the same J z value for any filling factor 
v, whereas for the short-range impurity, where the inter- 
valley scattering is significant, changing v may change 
(i) the number of transitions that are mixed and (ii) 
their J z values. This apparently affects the number of 
branches of bound states that are formed. Another dif- 
ference between the (5-function and Coulomb impurities 
is that there are degenerate states for the Coulomb im- 
purity but not for the (5-function. This is because the 
Coulomb impurity is sublattice symmetric, whereas the 
5-function impurity introduces an asymmetry between 
sublattices and hence valleys, thus breaking the SU(4) 
symmetry. Notice also, that the intervalley excitations 
K t> K' (f|"^-IJ-) are admixed by a short range potential 



to optically-active intravalley excitations opening a pos- 
sibility to introduce pseudospin flips by dipole photons. 

We believe that the application of the continuous de- 
scription to an extremely short-range <5-impurity system- 
atically underestimates the impurity effect. Indeed, only 
a few impurity interaction channels are open in a given 
LL in graphene as discussed above. Therefore, our pre- 
dictions for the energies of bound states are mostly of a 
qualitative nature. It would be interesting to compare 
these with results of the tight-binding approach incorpo- 
rating both the impurity and the Coulomb interactions; 
this may be a subject of a separate study. We wish to 
stress that, on the other hand, the symmetries of the ex- 
citonic states and the optical selection rules leading to 
dark and bright states established in this paper should 
remain valid in other schemes. We believe that these 
qualitative predictions can be tested well in experiments. 

A companion approach to studying the effect of short- 
range impurities would be the study of finite-range impu- 
rities, e.g., those which arise due to a disorder potential of 
Gaussian shape. Then the tuning of the height and width 
of the Gaussian would allow a continuous study of the 
transition from finite- to zero-range. Furthermore, one 
might then be able to quantitatively describe impurities 
induced by the presence of a suitably chosen substrate. 
We emphasize that the theoretical approach outlined in 
this work is equally applicable to such Gaussian disorder. 



Appendix A: Derivation of single particle 
Hamiltonian 



neighbor vectors as indicated in Fig. [T] and afffi) the op- 
erator which creates an electron on the A(B) sublattice. 
The Fourier transforms of the annihilation operators are 



«R, = —f= ^2 exp (ik ■ Rj) a (k) , 



(A2) 



J Rl+T„ 



-^=$^exp[tk-(R,+T TO )]&(k), (A3) 

k 

where N c is the number of unit cells in the crystal. Sub- 
stituting these into Eq. (|A1|) and expanding about the 
K, K' minima (the continuum approximation) yields 

™o = l E [(e i(K+q) ^<4 (q) & k (q) + H.c.) (A4) 



V A = 



W 
Nc. 



+ (V( K '+<^<4(q)frK'(q)+H.c. 

E [ a K (q) °k (q') + 4c (°0 »k' (q') 



(«k (q) ok' (q') + h. 



(A5) 



where, e.g., ok (q) = a (K + q). Since |q| <C 1/a, we can 
expand in (JA4I) to linear order in |q|. The summation 
over Tj is then performed using for graphene J^ ■ e 



iK-Tj 






iK -Tj 



iK'- 



and J2j e 3T i = 
tj = &a(i,—l). This gives 



Yaa*{i,l), 



The Hamiltonians, H^, for a single electron in the pres- 
ence of a (5-function impurity on a /i = A, B sublattice 
site at the origin, are derived below within the tight bind- 
ing and continuum approximations. The approach fol- 
lows that used in Ref. |2| to derive the free electron Hamil- 
tonian, Hq. The results are in agreement with the expres- 
sions obtained in Ref. [22| using a different approach. For 
an impurity at the origin and on the A sublattice, the 
tight binding Hamiltonian is 



H A = t 



Rj.Tj 



aLb Ri+Tj 



'R, 



a R; 



Wala 0l (Al) 



Ha = rio + Va- Here t is the nearest neighbors hopping 
energy, R^ the real space lattice vectors, Tj the nearest 



V3 



n = — ai^[(aq + ^(q) aK (q)+H.c.) (A6) 



q+ a K> (q) btc (q) + h.c. 



where q± = q x ± iq y and a — e 27 ™' 3 . To return to real 
space, we define the field operator ^ak. (r) according to 

a K (q) = Jdv e-( K +i' r * AK (r) (A7) 

and ^ak' (r) , f^K (r) > *bk' (r) similarly. Substituting 
these expressions into Eqs. (|A5[) and (IA6[) and converting 
the sum over q into an integral gives 



10 



Ho 



hvv, 



A N c . 



S_dq_ 

(2n) 2 

S_dq_ 
(2tt)* 



dv / dv' 



aq^ 



^(K+q). (-,-') ^ (r) % ak (r , } + jj_ a 



(A8) 



q + e 



(K'+q).(r-p')$t 



AK' ( r ) *BK' (r') + H.C. 

^ Jdrfdr' [(e^+^e-< K+ «> r '^ AK (r) * AK (r') + {K -j- K'} 

(K+q)r e^( K ' +q ') r '^t iK (r) * AK < (r') + {K «=» K'} 



Sdq' 

(2' 



(A9) 



where S is the area of the graphene sheet and v F = y/3at/2h, the Fermi velocity. The momentum operators are 

introduced by observing g±e ,q ( r ~ r ) = ^p±e* q '^ r ~ r > and replacing in the continuous approximation the canonical 
momentum p = —ihS7 by the kinematic momentum II — p + -A for the case of a magnetic field. The result may be 
expressed as 



H A 



/* 



/ *ak ( r ) \ T / V*(r) w F n_ VS(r) 



<K ( r ) 



V 






W/ 



w F n+ ooo 

F<5(r) F<5(r) v F U 
\ w F n_ 




*ak (r) 
a** BK (r) 

*AK' (r) 
-*bk' (r) 



(A10) 



which gives Eqs. (flj and (|2J). In Eq. (IA10p . V = \fiWa 2 /2 and the superscript T denotes transposition. Analogous 
derivation may be repeated for the case when the impurity is on the B sublattice chosen to be at the origin again, 
so the nearest neighbor vectors are given by the 8j indicated in Fig. [TJ whilst the valley definitions remain the same. 
This yields [cf. Eq. ©] 



Hi 



dv 



t 4w\ / 

«$BK ( r ) 

*U' (r) 





v F Il A 







« F n_ o 

VS(t) o 

o o 

-aV6(r) v F U. 




-a*V5{v) 

« F n + 

V«J(r) 



/ ®ak (r) 

a*$ BK (r) 

*AK' (r) 

V -*bk' (r) 



(All) 



r 



Note that different bases are used in Eqs. (|A10[) and 
(|A11|) due to the translation; they are related through 
*(r) = $(r-T i ). 



Appendix B: Higher order terms 

Here we estimate the relative magnitudes of higher- 
order energy corrections coming from the interaction 
with the 5-impurity. Suppose we have an electrically neu- 
tral system so that the lower cone is completely filled and 
the upper cone is completely empty. Let's denote a single 
electron state by \nsTm) = c}j Srm |0). As an example, we 
will suppose the impurity is on the B sublattice and cal- 
culate the ratio of the first and second order corrections 
to the single particle energy of an electron in the n = 
LL with, say, spin, pseudospin equal to J,, f|\ As indicated 
by the form of the impurity matrix element in Eq. © , an 
electron in the n = LL will only interact with the impu- 
rity if m = 0. Hence we consider the state 1 4-1*1" 0). The 
results are qualitatively the same for other scenarios. The 



The second order correction is 
E (2) _ y^ y^ i Kj3 oito 



\ V ni1tm i2 



IV. 



nlij.m<2 



n— 1 m— 



£CH1T ~ £r4iy 



(Bl) 



This represents virtual processes where an electron spon- 
taneously hops to a higher LL and then decays back 
to its original state. The sum over the LL index n is 
cut off at n c , which corresponds to the energy limit of 
where the linear dispersion relation ceases to hold. Using 

-ii nllYm c (— i) m V i -i i nlJLm c (—i) m V 



we obtain for Eq. (p3T 

E (2) 



*<4fro 



' 2V2ir£ 



V 2 



4V27r 2 hv F & 






(B2) 



B „=i 



This gives the ratio of the first and second order terms 

as 



E< 2 
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"c 



£W 2y/2Trhv F e B 4-| V™ 



E 
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-7=- ( B 3) 



first order correction is E^ = VboiIo = V/2tt£ 2 b ~ B. tion behaves as 



It can be shown that the ratio for the M order correc- 
r**s / \ M—x 

E (M) J y ^ ^ yM-Xg(M-X)/2 



p) 



(-x-Y 



11 



so that the perturbation theory is applicable when the 
dimensionless parameter V/hv^i-B <C 1. Note that in 
graphene this parameter increases with magnetic field B 
making the formal limit B — > oo singular; this is un- 
like the situation in the 2DEG, where the corresponding 
parameter is independent of B. In order to obtain a nu- 
merical value for Eq. (IB3|) . we need to estimate n c . Since 
the linear dispersion relation holds provided |q| <C 1/a, 
n c is defined according to hio c ^fn~ c = hvp/a. This yields 
n c sa 360 for B = 15 T. Using these parameters and tak- 
ing W = 1 eV for the impurity potential, we estimate 
\EW\/\EW\ rj 0.05. This shows that second and higher 
order terms in the perturbation expansion can be ignored 
for not too strong impurities and not in too strong fields 
B. 



Appendix C: Self-energy corrections 



f.V 



f.V 



5e^ 



f.v 



, > JlU! c 



w 



hu c + A% n 



-6- 



$tN &N = $W ~ Se^f 



FIG. 8: (Color online) Diagram illustrating the self-energy 
corrections, which contribute to the cyclotron resonance 
renormalization (increase). Note that "excitonic" corrections 
r to the cyclotron resonance energy due to the Coulomb ver- 
tex are not shown here. 



Kohn's theorem^ states that for electron systems with 
parabolic dispersion E(p) — p 2 /2m, the observed cy- 
clotron resonance energy is independent of the e-e in- 
teractions and is equal to the bare non-interacting value 
Huj c = heB/mc. The theorem does not hold in graphene 
with massless electrons. Formally, in conventional sys- 
tems with a constant charge-to-mass ratio, irrespective of 
dimensionality (2D or 3D), (i) the Coulomb vertex cor- 
rection T < to the energy [the e-h "exciton" lowering 
of energy analogous to Eq. ([ID])] is precisely compensated 
by (ii) the positive difference between the self-energy cor- 
rections A > to the final electron and initial hole states; 
the latter are due to the exchange with electrons in filled 
states. In graphene, the single particle energies also have 
self energy corrections, since each electron (hole) inter- 
acts via the exchange attractive (repulsive) interaction 
with all other electrons in the Dirac sea with same spin 
and pseudospin. However, this difference is not compen- 
sated by the excitonic effects in graphene, A + V > 0, 
leading to the Coulomb renormalization of the bare cy- 
clotron resonance energy huj c — > huj c = Hu> c + Shw c with 
5hoj c = A + T rj 0.7Eq, which may be absorbed in the 
renormalization of the Fermi velocity 13 ! 32 

Here we provide some details about the self- 
energy Coulomb correction to single-particle energies in 
graphene. The situation is represented schematically in 
Fig. [5] Such corrections are negative, of the order of Eq 
and decrease in magnitude as the LL number of the state 
in question increases. Specifically, the exchange interac- 
tion energy for an electron in state Afimi and a hole in 
state Mim-i is given by the matrix element U^ ™\m\ ml > 
as seen in Eq. (jlOp . Analogously, the interaction energy 
due to exchange between two electrons in states, Mm and 
M'm! is given by — U$£jj?$, where the sign change is 
due to switching h — > e. Hence the correction to the en- 
ergy of an electron in the LL n with spin s, pseudospin r 
and oscillator quantum number m due to exchange with 
electrons in the LL n', which has a fully filled same s,t 



sublevel, is 



Se(n,n') = -l^nrlm'n" 



TlTSm 
TSm' 



(CI) 



which may be written in terms of 2DEG matrix 
elements^ We remark that this expression is indepen- 
dent of the values of m, s and r. The complete self- 
energy correction for an electron in the LL n is given by 
S(-M = ~}2 n i 5e{n,n'), where the range over which n! is 
summed depends on the filling factor and the values of 
s, t. The renormalised single particle energies in Eq. Q 
are then given by ej^ — e^ + 5ej^. 

As a means of illustration, we shall give the self-energy 
corrections to the energies of those excitations mixed for 
the (A, <7 + , v = —1) system, shown in bold (red) in 
Fig. [5J Let Aj^f denote the self energy correction to the 
bare cyclotron resonance for an excitation with a hole in 
state Af and an electron in state A/"' '. The four excitations 
with no spin or pseudospin flip have the same corrections. 
Specifically for st = tJJ.,|f|\|JJ. 






A 0sr = 



A 0sT 

5(qst — Se-i s 



(C2) 
(C3) 



ST[Se(0,-n)-5e(-l,-n)] 



n=l 



e ^ r(n + |) y^ + 3 



8V^ 



n=l 



2n-l 



where T is the T-function and n c is the cutoff value for 
the LL index, as described in Appendix [B] One may 
show that A JY ~ lnn C) — so that for large enough values 
of n c , the result hardly changes for small alterations in 
n c . For B = 15 T, n c « 360 as before and A^J w 
1.43£'o'^ For the two excitations with a pseudospin flip, 
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Aq|| = A_°j{ = AjJJ + 0.25£ - This difference in self- 
energy corrections means that for excitations containing 
i. excitons with no spin/pseudospin flip and ii. pairs of 
excitons mixed by the direct interaction, the band width 
is greater by 0.25-Eo, than for excitations containing only 
one of these types of excitons. This leads to different 
band widths for different filling factors v — 0,2 and v — 
—1,1 (see beginning of Sec. lHI Al) . 
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A coherent superposition of such pair excitations consti- 
tutes collective excitations \N1N2Jz), see Eqs. (0) and (|18l) . 
This value differs slightly from the one we gave in Ref. [2U 
There the self energy correction due to exchange inter- 
actions with the lower cone only was given as 0.92_Bo- 
This corresponds to a total self energy correction, includ- 
ing exchange with both the lower and upper cones, of 
A(n c = 1870) = 0.75E +0.92E = 1.67£ . It is this value 
which should be compared with A(n c = 360) = 1.43-Eo 
quoted above. The remaining difference is because in this 
paper we use n c = 360, whereas in Ref. [2l] we took 
n c = 1870 as used in Ref. Ha . 



